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Abstract: Cosmic Microwave Background (CMB) polarization B-modes induced by Faraday Rotation (FR) can provide 
a distinctive signature of primordial magnetic fields because of their characteristic frequency dependence and because 
they are only weakly damped on small scales. FR also leads to mode-coupling correlations between the E and B type 
polarization, and between the temperature and the B-mode. These additional correlations can further help distinguish 
magnetic fields from other sources of B-modes. We review the FR induced CMB signatures and present the constraints on 
primordial magnetism that can be expected from upcoming CMB experiments. Our results suggest that FR of CMB will 
be a promising probe of primordial magnetic fields. 

PACS Nos.: 



1. Introduction 

Most of the ongoing theoretical and observational effort in 
cosmology is driven by its three most famous mysteries: Infla- 
tion, Dark Matter and Dark Energy. In comparison, the unex- 
plained origin of large scale cosmic magnetic fields may seem 
like many other problems in astrophysics - technically difficult 
but of lesser fundamental interest. This may, to some extent, 
be the case for magnetic fields in mature galaxies, where they 
could be generated through a dynamo mechanism [1]. How- 
ever, explaining their presence in young protogalaxies, clusters 
and, possibly, in the intergalactic space [2], is more challeng- 
ing. Evolution of magnetism in forming cosmic structures is 
highly non-linear which makes it difficult to conclusively rule 
out the existence of yet unknown astrophysical mechanisms 
that could generate magnetic fields on large scales and high 
redshifts. Another possibility is that magnetic fields existed 
before astrophysical structures began to form [3]. A primor- 
dial magnetic field (PMF) could be produced in the aftermath 
of cosmic phase transitions [4, 5] or in specially designed infla- 
tionary scenarios [6, 7]. Measurements of cosmic microwave 
background (CMB) temperature and polarization could deci- 
sively prove their primordial origin if they contained magnetic 
signatures present at the time of last scattering. A discovery of 
PMF would have profound implications for our understanding 
of the early universe, with critical insights into fundamental 
problems such as the matter-antimatter asymmetry [8]. 

A stochastic PMF influences CMB observables in several 
ways. Magnetic stress-energy perturbs the metric which leads 
to CMB anisotropics, while the Lorentz force deflects moving 
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electrons and protons coupled to photons. Recently, it has been 
suggested [9] that small-scale fields can appreciably alter the 
recombination history and, consequently, the distance to last 
scattering, because of the enhanced small scale baryonic inho- 
mogeneities. Here we focus on another signature of PMF - the 
Faraday Rotation (FR) of CMB polarization. 

FR produces a B mode type polarization with a character- 
istic spectrum [10, 11] as well as non-trivial 4-point correla- 
tions of the CMB temperature and polarization. In [12], we 
examined detectability of PMF using different correlators and 
evaluated their relative merits. We fond that a Planck-like ex- 
periment can detect scale-invariant PMF of nG strength using 
the FR diagnostic at 90GHz, while realistic future experiments 
at the same frequency can detect 10^^" G. This is comparable 
or better than other CMB probes of PMF, and using multiple 
frequency channels can further improve on these prospects. 

2. B-modes from Faraday Rotation 

At a given direction n on the sky, CMB is characterized by 
its intensity and two additional Stokes parameters, Q and U, 
quantifying its linear polarization. While Q(n) and U{n) are 
the quantities that experiments directly measure, their values 
depend on the choice of the coordinate axes. Instead, it has be- 
come customary to interpret polarization maps by separating 
them into parity-even and parity-odd patterns, or the so-called 
E and B modes [13, 14]. Existence of intensity fluctuations 
at last scattering implies generation of E modes, which by 
now have been observed and found to be consistent with the 
spectrum of temperature anisotropics. On the other hand, B 
modes would not be generated at last scattering unless there 
were gravitational waves or other sources of metric perturba- 
tions with parity-odd components such as cosmic defects [15] 
or magnetic fields [16]. Weak lensing (WL) of CMB photons 
by the large scale structures along the line of sight distorts 
polarization patterns generated at last scattering and converts 
some of the E mode into B modes, which is expected to be 
measured with upcoming CMB experiments. 

A primordial magnetic field present at and just after last scat- 
tering will Faraday-rotate the plane of polarization of the CMB 
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photons. The rotation angle along fi is given by 
3 



a(n) 



f B • dl 



(1) 



where f = nearCL is the differential optical depth, rig is the 
line of sight free electron density, <tt is the Thomson scat- 
tering cross-section, a is the scale factor, Ao is the observed 
wavelength of the radiation, B is the "comoving" magnetic 
field, and dl is the comoving length element along the photon 
trajectory. 

Statistically homogeneous, isotropic and Gaussian distributed 
stochastic magnetic fields can be characterized by a two-point 
correlation function in Fourier space 



{b,{k)b,{k')) = (27r)3,5(3)(k + k')[(% - k,k,)S{k) 



(2) 



where S{k) is the symmetric magnetic power spectrum, and 
where we omit the anti-symmetric contribution that quantifies 
the amount of magnetic helicity because only S{k) contributes 
to the CMB observables evaluated in this paper. The shape 
of S{k) depends on the mechanism responsible for production 
of PMF and generally can be taken to be a power law up to a 
certain dissipation scale: 



S{k) 



fc > fcdiss 



(3) 



The dissipation scale, fcdiss, should, in principle, be dependent 
on the amplitude and the shape of the magnetic fields spectrum. 
We assume that fcdiss is determined by damping into Alfven 
waves [20, 21] and can be related to Bctf as 



IMpc" 



10-'^Gauss\ 



(4) 



eft 



where B^s is defined as the effective homogeneous field strength 
that would have the same total magnetic energy density. B^s 
is related to the fraction of magnetic energy density to the total 
radiation density, ils-y, via [11] 



Boff = 3.25 X 10"^ V^Tb^ G auss . 



(5) 



The generation of CMB polarization and the PR happen con- 
currently during the epoch of last scattering. However, as we 
have shown in [11], assuming an instantaneous last scattering, 
i.e. that E modes were produce first and subsequently rotated 
by PMF, results in relatively minor inaccuracies. In this ap- 
proximation, the relation between the spherical expansion co- 
efficients of the E, B and a fields can be written as 



B, 



LAI linii 



lin 



kaLM 



(6) 



where A^/^^^ is defined in terms Wigner i-j symbols [17]. We 
note that B modes from WL can also be schematically written 
as (6) but with a different mixing matrix A^f|f,(^ . Importantly, 
the mixing matrix for WL has a parity opposite to that of PR 
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Fig. 1. The CMB B-mode spectrum from Faraday rotation in the 
case of a nearly scale-invariant magnetic spectrum with 2n = 0.1 
and Q.B-1 = 5 x 10^* (solid red), and a causal spectrum with 
2n = 5 and Sls^ = 10"^, at 61 GHz. The black short-dash 
line is the input E-mode spectrum, the black dash-dot line is the 
contribution from inflationary gravitational waves with r = 0.1, 
while the black long-dash line is the expected contribution from 
gravitational lensing by large scale structure. 

so that the two rotations are orthogonal to each other, making 
it possible to reconstruct them separately. 

In Pig. 1 we show the B-mode auto-correlation spectra due 
to PR by stochastic magnetic fields with two different spec- 
tra. One, with 2?! « 0, corresponds to nearly scale-invariant 
PMP generated via an inflationary mechanism [6, 7], while the 
other, with 2?! = 5, corresponds to PMP produced causally 
in phase transitions [22, 23], e.g. at the time of electroweak 
or QCD symmetry breaking. Also shown is the E mode auto- 
correlation spectrum which acts as a source for the PR B modes, 
as well as B modes from inflationary gravitational waves with 
r — 0.1, and the expected contribution from WL. 

The PR induced B mode spectra have certain characteris- 
tic features. In the case of the nearly scale invariant magnetic 
spectrum, the spectrum is oscillatory. In fact, the shape of the 
B-mode spectrum mimics that of the E-mode, except for the 
lack of exponential damping on small scales. This is because 
E modes are suppressed by the Silk damping, while PMP can 
remain correlated on small scales. The damping of the PR in- 
duced B-mode power is due to averaging over many random 
rotations along the line of sight. This translates into a 1/^ sup- 
pression of the angular spectrum, i.e. asymptotically we have 
p(jBB ^ /2n-i large I. In the case of 2n = 5, expected 
for causally generated magnetic fields [23], the PR produced 
B-mode can dominate the signal at high I. It is thus interesting 
to consider possibihty of future B mode experiments specially 
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Fig. 2. The minimum detectable magnetic field amplitude Qb-^ 
as a function of the magnetic spectral index 2n for the three 
estimators, BB (solid red), EB (dashed green) and TB (dotted 
blue). The top panel is for El, the middle panel is for E2 and the 
lower panel is for E3. 

designed to look for cosmological signals at sub-arcmin scales. 

3. Signal in mode-coupling correlations 

Spatially dependent FR produces additional non-Gaussian 
signatures in the CMB polarization. Namely, a particular real- 
ization of the FR distortion field correlates the respective Leg- 
endre coefficients i?/,„ and S/'m'- In fact, as shown in [17], 
it is possible to reconstruct the distortion field a{h) from spe- 
cially constructed linear combinations of products EimBi'm,. 
The additional correlations induced by FR also manifest them- 
selves as connected 4-point functions of the CMB, which, in 
turn, provide a measurement of the distortion spectrum 
[18, 19]. In principle, one can construct four estimators of 
the distortion spectrum, based on products of two CMB fields 
one of which contains polarization: TE, EE,TB, and EB. 
Of these four, the first two receive a large contribution to their 
variance from the usual scalar adiabatic Gaussian perturbations 
which makes it harder to find the FR signal. In [12], we stud- 
ied the last two, i.e we considered estimators based on 4-point 
correlations (EBEB) and {TBTB), and compared them to 
the 2-point function (BB) . We asked which of the estimators 
has the highest signal to noise for several types of magnetic 
field spectra and for a range of experimental sensitivities. We 
have fully accounted for the contamination by weak lensing, 
which contributes to the variance, but whose contribution to 
the 4-point correlations is orthogonal to that of FR. 

To forecast the detectability of FR we looked at three exper- 



imental setups: a Planck-like satellite [24] (El), a ground- or 
balloon-based experiment realistically achievable in the next 
decade (E2), and a future dedicated CMB polarization satellite 
(E3). The forecasts depend on the fraction of the sky covered 
by the experiment, /sky, which is close to unity for El and E3, 
and will be smaller for E2. We present our forecasted bounds 
on the PMF fraction in Fig. 2 subject to specifying /sky, which 
only appears under a quartic root in the bounds on Bcb- 

As Fig. 2 demonstrates, FR will be a very promising diag- 
nostic of PMF. In particular, future generation of sub-orbital or 
space-based CMB polarization experiments will be able to de- 
tect scale-invariant magnetic fields as weak as lO^^^G based 
on the measurement at 90 GHz frequency. Measurements at 
multiple frequencies can further significantly improve on these 
prospects. 

The relative strengths of the three estimators, demonstrated 
in Fig. 2 can be understood as follows. Generally, the EB 
and TB estimators have a larger number of independent modes 
contributing to the signal than the B-mode spectrum. Thus, in 
principle, it is not surprising if they result in a higher signal 
to noise. However, whether that is the case depends on the 
experimental noise level, and the distribution of power in the 
given combination of CMB fields and in the magnetic field. 
For a scale-invariant PMF spectrum, the i?-mode is essentially 
a copy of the E-mode, with most of the B-mode power be- 
ing on scales where the £'-modes are also most prominent. 
This results in a strong correlation between E and B for scale- 
invariant fields. In the case of the TB correlation, the under- 
lying T and E (B is obtained by a scale-invariant rotation of 
E) fields peak on rather different scales. Namely, T peaks at 
i 200 while E peaks at £ ~ 1000. In other words, the 
intrinsic correlation between T and E is already suboptimal, 
translating into a lesser correlation between T and B. Thus, 
for experiments with sufficiently low noise levels, such as El, 
E2 and E3 considered in this paper, the EB estimator performs 
better than TB for scale-invariant fields. This would not nec- 
essarily remain true if polarization measurements had a signif- 
icantly higher experimental noise. 

For causally generated PMF with steeply rising ("blue") spec- 
tra, as in the 2n = 5 case, the FR power is concentrated on very 
small scales, far away from the scales at which any of the un- 
rotated CMB fields have significant power. This means that 
the B-modes in the observable range are obtained either by 
a rotation of i?-modes far away from their peak power scale, 
or by a rotation of peak i?-mode by a negligible angle. This 
means that E and B fields peak at very different scales, with 
their correlation being close to zero over the observable scales. 
In this case, we see that the S-mode spectrum, i.e. the BB 
correlation, has the highest signal to noise. 



4. Discussion and Outlook 

When interpreting the forecasted bounds on the magnetic 
field energy fraction or the effective magnetic field strength in 
Fig. 2, several points must be kept in mind. First, the con- 
straints are obtained after setting the dissipation scale to be 
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given by Eq. (4). For scale-invariant fields it makes no dif- 
ference, as in this case most of the signal is on scales larger 
than the magnetic dissipation scale, and fcdiss does not con- 
tribute to the normalization of the spectra when 2n 0. In 
fact, for scale-invariant fields, the effective field Bcs defined 
via Eq. (5) is the same as the commonly used B\, which is 
the field smoothed on a given scale A. Thus, our forecasts of 
the minimum detectable BcS for scale-invariant fields can be 
directly compared to most other bounds in the literature. 

CMB is less sensitive to magnetic fields with blue spec- 
tra because most of the anisotropies are concentrated on very 
small scales. This is what Fig. 2 is showing too. Here we note 
that Fig. 2 assumes that the spectrum will keep rising at the 
same steep rate (2?! = 5) all the way to the dissipation scale. 
On the other hand, simulations [23] suggest that the spectrum 
becomes less steep, with 2n' = 3 over some intermediate 
range kj < k < fcdiss, implying a smaller net magnetic en- 
ergy fraction flBf Big Bang Nucleosynthesis constrains this 
fraction to be less than 10% [27-30]. Note that the expected 
bound from Planck (El) in Fig. 2 for 2n = 5 will be at least an 
order of magnitude stronger, while E2 and E3 will improve on 
the BBN bound by two orders of magnitude. Here it is worth 
keeping in mind that our bounds are on the magnetic field con- 
tribution at the time of last scattering. While it is expected 
that the fields remain effectively frozen-in between the time of 
nucleosynthesis and last scattering, with a relatively slow time 
evolution of the dissipation scale, this is still an approximation. 
In any case, it is interesting to know how good the resolution of 
future B mode experiments can be, since the PR contribution 
from causally generated PMF appears to dominate over other 
cosmological sources on small scales. 

In the case of scale-invariant fields, existing bounds on the 
magnetic field strength from WMAP are at a level of a few 
nG [31]. These bounds are based on the anisotropies induced 
by the metric fluctuations sourced by magnetic fields, and ig- 
nore the FR effect. In Refs. [11, 32] the WMAP bound using 
FR was obtained at the 10~^ G level. As one can see from 
Fig. 2, Planck (El) can almost match today's bounds for scale 
invariant (n = 0) fields using the EB estimator at only one fre- 
quency, while future probes, such as E2 and E3, can improve 
the bounds by an order of magnitude! This suggests that the 
mode coupling estimators of FR can be a very powerful direct 
probe of scale-invariant PMF. 

Finally, as mentioned already, the bounds are based on us- 
ing a single frequency band, while using several bands can fur- 
ther improve the constraints. While our estimates look quite 
promising, they are still preliminary and ignore the potentially 
devastating foreground effects. Prior to reaching our detectors, 
CMB must pass through the magnetic field of our own galaxy 
and will experience FR in which B modes are produced. It re- 
mains to be shown to what extent the FR due to the galaxy can 
be subtracted from the cosmological FR signal. 
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